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The interaction of whispering gallery modes (WGM) of optical microresonators with subwave- 
length imperfections has been studied both experimentally and theoretically. This interaction is 
responsible for the formation of spectral doublets in place of single resonance peaks, and for degrad- 
ing of Q-factors of the resonances. Within the currently accepted framework the spectral doublets 
OO ' are explained as a result of degeneracy removal of clockwise and counterclockwise WGMs due to 

' their coupling caused by defect-induced backscattering, while the degrading of the Q-factor is de- 

, scribed phenomenologically as an additional contribution to the overall decay rate of WGM due to 

coupling between WGM and radiative modes. Here we show that the existing understanding of this 
^ ' phenomenon is conceptually wrong and develop an exact theory of WGM interaction with a single 

1^ , defect, which provides a unified treatment for both aspects of this interaction explaining existing 

■ experiments and predicting new phenomena. 

m ■ 

^SJ \ Elastic (with no change in frequency) scattering of hght due to small (compared to wavelength) particles is one of 
the most fundamental and intensively studied optical phenomena. Its modern history began almost one hundred fifty 
years ago with the explanation of the blue color of sky in a series of papers by Lord Rayleighfl'l, where the now famous 
O i l/A'* cross section law, where A is the wavelength of light in vacuum, was derived. Since then it has been customary to 
^ ■ refer to processes of elastic interaction of light with subwavelength particles as Rayleigh scattering. Besides providing 
us with beautiful blue skies and red sunsets, Rayleigh scattering is important for a large number of fundamental 
^ optical phenomena as well as for numerous applications. Recent developments in optics and photonics have created 
^ \ new situations in which the manifestations of Rayleigh scattering are significantly modified. Particularly drastic 
^ . modification of this process is expected when light is confined in all three dimensions inside optical microresonators 
^ in the form of whispering gallery modes (WGM) [2] . Given the fundamental nature of this process it is not surprising 
^I^i; that it has attracted a significant amount of attention in recent years [1, [3, 0| • 

While whispering gallery modes can occur in various types of geometries [8| we will focus on spherical microres- 
onators. WGMs in this case correspond to Mie resonances @ with ultra narrow widths, 7/s ^ ujis, where ujis is the 
frequency of the mode, and respectively high (up to 10^ for silica microspheres fs^) Q-factors defined as Qis = ujis/^is- 
WGMs are characterized by polar and azimuthal indexes, I and m, and a radial number s determining, respectively, 
the angular and radial dependence of the fields in a spherical coordinate system centered at the sphere. The resonance 
' frequency w/s does not depend on the azimuthal number, which reficcts the degeneracy of the resonances due to full 
: spherical symmetry of the problem. WGMs are also characterized by their the mode volume, which can be very 
■rj" ■ different for modes with the same I but different m. Modes with the smallest volume correspond to \m\ = I, and s = 1 
in which case the field is concentrated mostly in the equatorial plane and at the surface of the sphere. Such modes 
are called fundamental (FM) and their interaction with defects is of the primary interest. 
00 This interaction causes two observable effects: (i) formation of spectral doublets in place of a single peak, and (ii) 

' reduction of Q-factors of WGMs below theoretically predicted limits 0, [l0| . In the existing approaches these two 
^ ' effects are considered separately, even though they are two manifestations of the same phenomenon. The accepted 
explanation of the spectral doublets is based on the hypothesis likely first suggested by D.S. Weiss et al. in the 
following form: "We have observed that very high-Q Mie resonances in silica microspheres are split into doublets. 
This splitting is attributed to internal backscattering that couples the two degenerate whispering-gallery modes 
propagating in opposite directions along the sphere equator" Q . This idea of the defect-induced backscattering was 
further developed in subsequent publications 4, 5, 7] and was claimed to be experimentally confirmed in Ref. 
In a recent paper by A. Mazzei, et al. [ll| double peak features in the spectra of microspheres were studied under 
conditions of controlled scattering, where the role of the defect was played by a tip of a near-field optical microscope. 
This work directly confirmed a connection between the interaction of WGMs with a single defect and the formation 
of the spectral doublets. 

The influence of defects on Q-factors of microresonators is usually prescribed to defect-induced coupling between 
WGMs and radiative modes and is taken into consideration phenomenologically by adding a "scattering" loss rate to 
the total losses of the resonator. In the case of disk resonators this rate was calculated in Ref. ^ under the assumption 
that the surface roughness couples WGMs with free space electromagnetic radiation. 
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FIG. 1: The coordinate systems used in calculations. The curved arrow schematically represents the counterclockwise funda- 
mental mode whose interaction with the defect (blackened small sphere) is under consideration. 



I. RESULTS 



In this work we develop an ab initio theory of interaction between WGMs of microspheres and a single subwave- 
length defect, which treats both aspects of this interaction in a unified way. We show that the generally accepted 
picture of interaction between FM and defects is conceptually wrong: the interaction cannot be described in terms of 
" backscattering" , and it does not result in coupling between counterpropagating " clockwise" (cw) and " counterclock- 
wise" (ccw) FMs. The theory predicts that the observed spectral doublets are actually a part of a triplet of peaks, 
the third component of which has not yet been found. It also provides an accurate "from first principles" description 
of the broadening of the resonances. Similar to the Rayleigh scattering of propagating waves, the solution of the 
single-defect problem presented here constitutes the first fundamental step toward a theory of interaction between 
WGMs and multiple defects. When this interaction is small, which is usually the case, the single-defect solution 
can be directly used to make conclusions about the role of multiple scatterers. Since the results of Ref. [llj directly 
confirm that the spectral features formed due to a single discrete scatterer are similar to those caused by distributed 
surface or volume disorder, the results presented here are relevant not only for interaction of WGMs with discrete 
non-interacting defects, but can also be used to understand effects due to continuously distributed nonuniformities. 

The main assumption of the theory is that the defect is small enough to be treated as a dipole. In this case the 
shape of the defect is irrelevant, and can be taken to be spherical. In this way the problem is reduced to describing 
two electromagnetically coupled spheres of radii Rg and Rd ^ Rq with refractive indexes n and Ud respectively, whose 
centers are positioned at a distance d from each other (see Fig. [T]). We will also assume that the defect lies in the 
plane of the FM in order to maximize the strength of the interaction, although the theory developed here can readily 
be generalized for arbitrarily positioned defects. The goal is to find electromagnetic field induced by this system in 
the presence of an incident wave (imitating a mode of a tapered fiber) which, in the absence of the defect, would 
have excited a FM with given polar number L. The incident, Ei„c, and induced, E^, fields can be presented as linear 
combinations of vector spherical harmonics (VSH) of the form: 



Einc = ^ VL,7n'NL,7n{r - ri) (1) 
m— — L 
2 oo / 

E« = EE E hlN,™(r-r,0 + 6«M,,„(r-r,)] (2) 

2—1 l — l ■m— — l 



where index i enumerates the spheres (i = 2 refers to the defect), ri is a position vector of the center oi i ~th sphere, 
M;_,„ and N;.,„ are the VSH of TE and TM polarization respectively as defined in Ref. [l^], and r?L,m are coefficients 
describing the TM incident wave of frequency uj, which in the coordinate system XY Z defined in Fig. [T]and used in 
all subsequent calculations have the following form 



,( lWL+^)H)^^/ (2L)! \1 cwFM 

^^•""^ ' 2^ ^ {L + m)\{L-my: \0 ccwFM 

(See details in Section ITTTl ) In order to find the induced field one needs to determine expansion coefficients a[*^ and 
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FIG. 2: Radiated energy of the microsphere-defect system with varying distance parameter S — {d — R{) — Rd)/Ro- 



b'fl^. For simplicity wc disregard the defect-induced coupling between TE and TM modes, and since the incident field 

is of TM polarization we set the TE coefficients — and solve for coefficients aj^^''. The dipole approximation 

for the defect is introduced by setting ^ = for all / > 1. The resulting system of equations for the scattering 
coefficients is solved exactly to yield: 



,(2) 



= < 



a 



(1) 

O^l VL,mC>LL 



(2)u,(l)j-l^ 



L;\m\<l (a) 

l^L;\m\<l (b) 
\m\ > 1 (c) 



(4) 



(5) 



where a*-^'^^ are the single sphere scattering parameters for the main sphere and the defect respectively, defined in 
Eq. \TE\ and \W\ of section IIIII The scattering parameter of the main sphere has poles at the complex- valued frequency 
of WGMs; we assume that the frequency uj of the incident field is in the vicinity of the pole w^"'' 



iri°'' corresponding 



to the frequency of our FM. Since the defect is assumed to be small so that UdUjRd/c ^ 1, its scattering parameter 
does not have any poles of its own. 

Parameters ^z"'™ in Eq. |4] and [5] describe the electromagnetic interaction between the spheres and are called 
translation coefficients. They appear when a VSH defined in one coordinate system needs to be expressed in terms 
of VSH defined in a system with a shifted origin [l^, [ij, [l^ . These coefficients depend on the translation vector 
ri — Y2 and the coordinate system used to define VSHs. In the coordinate system XYZ they are diagonal in terms 
of azimuthal number m, which reflects the fact that the polar axis Z runs along the line connecting the centers of 
the spheres, thus preserving the axial symmetry of the two-sphere structure. Eqs. 2] and [5] contain all the information 
about the electromagnetic field of the sphere-defect system. First of all, Eq. \^ shows that components of the initial 
FM with |m| > 1 are not effected by the defect. Formally, this result is a consequence of the translation coefficients 
^r'm being equal to zero when either m or /i exceeds either v or / fl5|. Physically, this result refiects the simple fact 
that a dipole can only produce a field with I — 1 and m = 0, ±1. Since in the XY Z coordinate system m remains a 
conserving quantity even in the presence of a defect there will be no coupling between the field of defect and WGMs 
with \m\ > 1. Thus, the expansion coefficients a|^2i with I = L and \m\ > 1 will produce a resonance at the original 



single-sphere frequency, which we will characterize by size parameter x]- 

„(o) 



The width of this 



resonance, which can be described by a dimensionless parameter j'"^'' — vf^ Ri)/c is also not affected by the defect. 
Eq|4^, on the other hand, shows that |m| < 1 components of the FM do interact with the defect, and that this 
interaction results in appearance of new resonance frequencies determined by poles of this expression. From the sum 
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over V we singled out a term with i/ = L, which in the frequency range around Lof^ gives the biggest contribution to 
the shift of the new poles from their single-sphere value. Discarding the rest of the sum (the resonance approximation) 
we can obtain an analytical expression for the positions of the new poles 



(0) 



(2L + l)Rld 



where 



^ r (0)i c 2fL,ra(k'-^^d) 



(6) 



3 (2L+l)i?5d^ J 



(7) 



is the standard polarizability of a small dielectric sphere, and function /l,™ {kf^dj for m — 0^ ±1 is defined as 
fL,rn{kd) = 



i-iy 



(i + l)(i + m2) 



1 + 



gL-i{kd) + \jL{L + 1)(1 - m2) + L'^—gL^,{kd) 



(8) 



with 



gL{kd) = 



(atanhp— p) 



where we used an asymptotic form of the Hankel function valid for / ^ kd. Eq. [6] predicts two new resonances, in 
addition to the original single sphere resonance, one for m — and another for m = ±1, both red shifted with respect 
to the initial frequency. Function fLm{kd) specifies the dependence of the resonance frequencies on m and distance 
d. The later is determined by the exponential decay of the spherical Hankel functions outside of the main sphere, 
which reflects the evanescent nature of the interaction with the defect. Thus, our theory predicts the existence of a 
triplet of peaks rather than the doublet expected in the current cw-ccw coupling picture. To verify this result we 
carried out numerical calculations of the frequency dependence of the energy emitted by the main sphere given by 
the standard expression + l)|a!,mP using complete Eq. [H For these calculations we choose L = 39 and 

take into account enough coefficients a/_m and terms in the sum over u to ensure convergence of the procedure, which 
was achieved with 1 < ^ < 50 and < 50. The results of these calculations are shown in Fig. [2] for different distances 
between the defect and the sphere so that one can see how the peaks shift toward the single-sphere resonance with 
increasing d and eventually merge with it. There are indeed three peaks, which are not seen at curves 1 and 2 because 
the third peak on these curves is out of the range of the figure. The to = resonance is shifted further from x^^^ and 
is weaker than the to = ±1 resonance, making it more difficult for experimental identification. We suggest, therefore, 
that the experimentally observed spectral doublets correspond to the original single-sphere resonance and the m — ±1 
resonance introduced by the defect. 



The validity of the resonance approximation depends on the convergence of the sum J^u^Li"^)'^'^ 
appearing in the denominator of Eq. [4^. In the limit v —^ oo we find 



(-l)''+ia(i)A-:A 



l.m 



{kdf 



Ro 
~d 



2u+l 



(9) 



which, given that Ro/d < 1, proves the convergency of the sum. However, in the case of small defects positioned close 
to the surface of the sphere, Ro/d differs from unity by a small amount and the convergence of the sum is slow. In this 
case the terms with p ^ L become important and should be taken into account. It can be shown that incorporating 
these terms does not change the form of Eq. [B]but renormalizcs the polarizability, which becomes: 



P = P 



1 



1 



TO 

T 



dRliRl 



1 



{kdf{d^~RlY 



+ 2 



(10) 



The renormalized polarizability acquires dependence on the distance d between the defect and the sphere, thereby 
affecting the relation between the frequency shift and the broadening of the defect-induced resonances which is revealed 
through variation of d. Indeed, in the absence of the renormalization both these quantities decrease with d by the 
same factor determined by the function fL.m, resulting in 5xL,m oc 5^L,m- The renormalized polarizability, however. 
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FIG. 3: Relative broadening versus frequency shift of the |m| = 1 resonance revealed through varying distance d. The points 
are obtained from the spectra computed at different d, while the line represents a fit with a quadratic polynomial 
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FIG. 4: Internal field intensity of the microsphere in the YZ plane at the frequency of the standard Mie resonance (left) and 
the defect induced resonance (right). 

also changes with distance. Since the frequency shift of the resonance is linear in p while the broadening is quadratic^ 
this effect must result in deviations from this linear dependence. In order to confirm this conclusion we used numerical 
spectra obtained for different distances in order to plot Sx^^m versus S^L.m for the \m\ = 1 resonance. The obtained 
data shown in Fig. [3] are found to be better fit by a quadratic rather than a linear function confirming this conclusion. 

With scattering coefficients known we can also compute the internal field inside the main sphere. Fig. [4] shows 
the variation of the field in the YZ plane of the XYZ system (which corresponds to the plane of the FM) obtained 
by varying radial and polar coordinates at the azimuthal angles (j) = Tr/2 and (f> — 37r/2. The computed field profile 
for the defect-induced peak demonstrates 2L oscillations and a drastic increase in intensity in the vicinity of the 
defect. At the frequency of a single-sphere resonance the situation is reversed: 2L oscillations, which are phase shifted 
compared to the defect-induced resonance are accompanied by a significant decrease in the field's intensity in the 
defect's proximity. 

Finally Eq. describes coupling of the FM to other WGMs, most important of which are terms with I < L. 
There are two reasons for this. First, modes with lower I and higher radial numbers can spectrally overlap with the 
Z = L, s = 1 mode [l^ and thus have a large effect on the field distribution. Second, these modes usually have 
lower Q-factors and contribute more significantly to the radiation losses of the system. Note that m-components with 
TO > 1, which are responsible for the resonance at the single-sphere frequency do not couple to any other modes, so 
that this resonance is not affected by the coupling to the low-Q WGMs. The effects of coupling to these modes at 
the frequency of the defect-induced peak is shown in Fig. [5l where we plot the spectrum of the radiated energy in its 
vicinity with and without contributions of terms with I ^ L. 
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FIG. 5: Radiated energy with (the curve with the higher peak) and without contributions from I ^ L. 

II. DISCUSSION 

The theory presented in the paper gives a complete picture of the interaction between WGMs and a single defect 
based on fundamental principles with no ad hoc assumptions, and replaces the currently accepted paradigm, which 
is proved to be inadequate. The results obtained on the basis of this theory show that a number of experiments 
previously " explained" within the cw-ccw coupling picture such as the backscattering experiment of Ref . \u\ must be 
reinterpreted. They also provide a natural explanation to other experimental results that the current paradigm was 
not able to explain in addition to predicting new effects that await experimental confirmation. 

To begin with, the developed theory gives a natural explanation of an asymmetry between the two peaks of a 
doublet, which was seen in all experimental observations of this effect, but most clearly in Ref. [ll|. Since, according 
to our calculations, the higher frequency component of the doublet corresponds to a single sphere resonance unaffected 
by the defect, this peak is supposed to be narrower than its counterpart and not to shift with the change in the position 
of the defect. This behavior is in complete agreement with observations of Ref. [11]. In addition, our theory predicts 
the existence of the third peak, which is, however, weaker than the other two making its experimental observation 
more difficult. It should be noted, however, that since within the prevailing paradigm the presence of the third peak 
was not expected, it is possible that more careful experimental observations will reveal its presence. Eq.[6]also predicts 
quite specific dependence of the peak's position and its broadening on the position of the defect and the polar number 
L and frequency w^""* , which also can be verified experimentally. 

We also found that while the main contribution to the width of the resonance comes from coupling to the dipole 
field of the defect, there is also an additional contribution from coupling to lower Q WGMs of the main sphere. This 
contribution is emphasized in Fig. [51 where we compare defect-induced resonance with and without terms with / ^ L. 
One can see that these terms make the resonance wider while increasing its height. This effect, which cannot be 
described by simply adding an additional loss term to the Q-factor, is relatively small for a single defect case, but can 
be expected to become more significant with multiple defects. 

Another important effect predicted in our theory is the position dependent renormalization of the polarizability of 
the defect. This effect results in a deviation from the linear dependence between the frequency shift of the defect- 
induced resonance and its broadening, which was recently observed in Ref. [ll| . The authors of that paper suggested 
that the position dependence of the polarizability could account for this finding, but were unable to explain its origin. 
In our theory this dependence appears naturally as a result of coupling between the defect and WGMs with high 
polar numbers I. 

An important characteristic of the interaction between the FM and the defect is the resulting distribution of the 
electromagnetic field along the surface of the sphere. The backscattering paradigm predicts the formation of standing 
waves with 2L oscillations of the field's intensity along the circumference of the FM. These waves are assumed to 
be due to interference of cw and ccw modes and are described by either sin- or cos-like behavior, depending on 
which component of the doublet is considered. Using Eq. [3] we can see that in XYZ coordinate system used in our 
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calculations these standing waves should be described as 



'SW — 




(L + m)!(L-m)! 



{2L)l 



■NL,™(r-ri) 



(11) 



where for the symmetric combination of the cw and ccw modes m takes on only even values for even L, and odd values 
for odd L; for antisymmetric combination the situation is reversed. Our results show that for both the defect-induced 
and single sphere resonance there are indeed 2L oscillations, which, however, do not have the form prescribed by 
Eq. [TlJ The field distribution at the frequency of the defect-induced peak is explained by the fact that the field at 
this frequency is mainly comprised of the m-components with |m| = 1. The field of these WGMs is characterized by 
L — \m\ + 1 = L oscillations for 9 changing between 0, tt giving their total number equal to 2L. These modes are also 
characterized by the enhancement of the field in the vicinity of 6* = 0, which explains a drastic rise in the intensity 
around the location of the defect. The field distribution at the single sphere resonance can be understood by noting 
that this field is comprised of modes with |m| > 1, which when added to the remaining \m\ < 1 components, would 
have produced a flat distribution of the intensity. Therefore, removal of these components obviously results in the 
decrease of the field around the defect and phase shifted oscillations elsewhere. The presence of these oscillations 
of the field's intensity demonstrate that one can explain experimental results of Ref. @ without reliance on the 
"backscattering" paradigm. 

Finally, we should comment on the relation between our results and the multi-defect problem. In the approximation 
of non-interacting defects, the field in the presence of multiple defects can be found as simple sum of fields due to 
each defect separately. Therefore, the generalization of these results would include finding the field distribution for 
a generic position of the defect relative to the plane of the FM. Since interaction with each defect can be considered 
independently we can always analyze it in a coordinate system with polar axis passing through the centers of the 
defect and the main sphere. Then, we can repeat all our calculations with only one adjustment: Eq. [3] needs to be 
generalized to incorporate an arbitrary inclination of the plane of the FM with respect to the coordinate axes. It 
is clear, therefore, that while the resulting field distribution and the heights of the resonance peaks will be different 
from the ones obtained here, the real and imaginary parts of the resonance frequencies will remain the same as in 
the single-defect case as long as all defects are identical [13]. It is also important to note that these quantities do not 
depend on the type of the initial FM, cw or ccw, which means that even in the presence of a defect or multiple defects 
there are still two modes originating from cw or ccw initial FMs, which remain degenerate. This conclusion has a 
number of far reaching implications and dispels another important myth of the cw-ccw coupling paradigm, which 
suggests that in the presence of the defects each resonant peak corresponds to a single non-degenerate mode. 



Results presented in this paper were obtained with the help of a number of qualitative and quantitative methods. 
We will start with the approach based on symmetry arguments, which provides a qualitative explanation of many of 
the findings. 



The logic behind the ccw-cw coupling paradigm is based upon an implicit assumption that the degeneracy between 
cw and ccw modes, which persists even in the absence of the complete spherical symmetry, is due to the rotational 
symmetry with respect to the axis perpendicular to the plane of the FM. Indeed, if this were the case, then any 
defect would have violated this symmetry, lifting the cw-ccw degeneracy and resulting in the spectral doublet. This 
assumption, however, is not correct, which can be immediately seen if one recalls that the group of rotations about a 
single axis is Abelian and, therefore, can only have one-dimensional representations. This means that axial symmetry 
alone cannot explain the cw-ccw degeneracy and one needs to invoke an additional symmetry such as inversion with 
respect to the azimuthal angle (j). In the case of a single sphere the inversion symmetry is "hidden" behind the more 
powerful spherical symmetry, but in the system with a defect it starts playing a significant role. Indeed, two interacting 
spheres, described in X'Y'Z' coordinates (Fig.[T]), which is not consistent even with the remaining axial symmetry of 
the system, still exhibits a symmetry with respect to replacement <j) ^ —0 if the X-axis of that coordinate system is 
chosen along the line connecting the centers of the spheres. Since this is the symmetry ultimately responsible for the 
cw-ccw degeneracy and it is not destroyed even in the system with the defect, the alleged coupling between cw and 
ccw modes cannot take place. 



III. METHODS 



A. Symmetry considerations and choice of the system of coordinates 
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The X'Y' Z' coordinate system is convenient to describe the FM excited in the main sphere shown in Fig. [T] In 
this case the respective field can be presented as a single vector spherical harmonic of TE or TM polarization [13] , 
meaning that the expansion coefficients of Eq. [1] takes a simple form r\L^rn — ^L.m instead of those given in Eq. [3l 
However, in the two-sphere problem this coordinate system is not consistent with the symmetry of the configuration, 
therefore, it is more convenient, following Ref. [13) [III > switch to a coordinate system with polar axis directed al ong 
the line connecting the centers of the spheres (designated as XYZ in Fig[T]). It is important to realize, however, [l^l 
that the field of this FM cannot be presented as a single VSH in the spherical coordinates based on the coordinate 
system XY Z . To obtain such representation we notice that this system is obtained from X'Y' Z' system by means 
of a rotation characterized by Euler angles a = tt/^^B = 7r/2, 7 = 0, where we are following notations from Ref. [isj . 
Now using the transformation properties of VSH [1^ we obtain the representation of the FM in the XY Z-h&sed 
spherical coordinates given in Eqs.[l]and [31 

The XYZ coordinate system reflects the presence of the axial symmetry of our configuration with respect to 
rotation about the axis connecting the centers of the spheres. Because of this symmetry, even though the modes of 
the sphere with the defect can no longer be classified according to the polar number, I, they still can be characterized 
by azimuthal number, m. Respectively, each of the m-components comprising the FM interacts with the defect 
independently making the analysis of the interaction simpler. 



B. Multi-sphere Mie theory 

In order to find the expansion coefficients of the induced field introduced in Eq. [2]we use the standard multi-sphere 
Mie theory In this approach, the field outside of the spheres is separated into incident field given by Eq.[T] 

and the induced field given by Eq|21 In addition the field inside the spheres is also presented as a linear combination 
of VSHs centered at each sphere 



[q« N,,™(r - r.) -f 4iM,,™(r - r.)] . (12 



In order to apply Maxwell boundary conditions on the surface on a sphere i the VSHs centered at different spheres 
must be rewritten in the coordinate system translated to the center of the i-th sphere. This is accomplished with 
the help of the addition theorem for the vector spherical harmonics [13, [13] , which introduces translation coefficients 
■^i m i^'^j ~ ^i) ^'^d i^^^j ^^^^1 which allows one to derive a system of equations relating expansion 

coefficients aj*^ and to the coefficients of the incident field jy^*^: 



= { vl% + E E - r.) + I^IX;::' (x, r, - r.)] (13) 

j^i I' ,in' j 

C = Cl^' E E - r.) + a^l,Bl:^' (x, r, - r.)] (14) 

j^i I' ,m' 

(i) (i) . 

where and Q are single sphere Mie scattering parameters for TM and TE polarizations respectively. For the 
TM polarization this parameter, defined in terms of dimensionless frequency parameter x — Rqlu/c, is given by the 
well-known expression [Tsj 

^(1) _ jiix)j;[xji{nx)] - n^ji{nx)-^[xji{x)] ^^^^ 
' ^' (^) ^ i^ji i^^)] ~ "^^ji ("2;) ■£ [xhi {x)] 

where ji{x) and hi{x) are Bessel and Hankel functions respectively. For the defect we only need I = 1 and if n^Xd ^ 1, 
where Xd = xRd/Ro the scattering parameter a\ does not have any poles and can be approximated as 



„(2) 



Explicit expressions for translational coefficients Ai '™ (vj — r.;) and _B; ^ {vj — r.;), which describe optical coupling 
between the spheres via modes of the same or different polarizations respectively, can be found, for instance in 
Ref. [l^ [13, [la]- Important property of the translation coefficients is that they take a diagonal form in m if the 
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translation vector is parallel to the polar axis of the coordinate system used to define spherical coordinates. This 
significantly simplifies the equations for expansion coefficients eliminating summation over the azimuthal number and 
decoupling equations for coefficients with different to. We take advantage of this property by working in the XY Z 
coordinate system of Fig. [T] 

For WGM with I ^ 1 the cross-polarization translation coefficients are usually much smaller than their same- 
polarization counterparts. Since we assumed incident wave to be of TM polarization we can set bi^m = in Eq. 1131 
obtaining as a result a closed system of equation for coefficients ai^m- 

C. Dipole approximation 

The field of the dipole is described by VSHs with 1 = 1. Therefore, we introduce the dipole approximation by 
assuming that ap2i = for Z > 1. This reduces the system of Eq. [T3]to a simpler form 

C = «PHe + «S„4'r(ri-r,)} (17) 
ail = (ri - r2) (18) 

which can be solved exactly by multiplying Eq.[T7]by (— 1)^+'A'^'™ and summing over I. Substituting Eq.[TH]into the 

resulting expression, we obtain a closed equation for the quantity X]i/(^l)^^'^al.lA!j''™, which can be easily solved. 
As a result we arrive at Eq. 2] and [5] of Section [D 
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